Time-symmetric fluctuations in nonequilibrium systems 
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For nonequilibrium steady states, we identify observables whose fluctuations satisfy a general 
symmetry and for which a new reciprocity relation can be shown. Unlike the situation in recently 
discussed fluctuation theorems, these observables are time-reversal symmetric. That is essential for 
exploiting the fluctuation symmetry beyond linear response theory. Besides time-reversal, a crucial 
role is played by the reversal of the driving fields, that further resolves the space-time action. In 
particular, the time-symmetric part in the space-time action determines second order effects of the 
nonequilibrium driving. 
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Nonequilibrium statistical mechanics is still very much 
under construction. Its very rich phenomenology reach- 
ing from some of the great unsolved problems of classical 
physics such as turbulence, to the interdisciplinary helds 
of biophysics and complex systems, has not yet found a 
sufficiently general and powerful theoretical foundation. 
Compared with equilibrium statistical mechanics, we lack 
a global fluctuation theory of far from equilibrium sys- 
tems, we do not have a controlled perturbation theory 
and even many conceptual debates have not found their 
generally accepted and logical conclusions. In contrast, 
the Gibbs formalism of equilibrium statistical mechanics 
yields general identities and inequalities connecting ther- 
modynamic quantities and correlation functions. They 
relate the response of a system to an external action with 
the fluctuations in the system of the corresponding vari- 
able. That is partially inherited by irreversible thermo- 
dynamics, as in the fluctuation-dissipation theorem, 
that works close to equilibrium but is mostly limited to 
the so called linear response regime; away from equilib- 
rium no such general relations exist. 

In recent years however 0, 0, El EJ 0> 13 a broader 
range of thermodynamic and fluctuation relations have 
become available also further away from equilibrium and 
there has been a general sentiment that these are impor- 
tant steps in the construction of nonequilibrium statis- 
tical mechanics. In what follows we limit ourselves to 
the steady state. A fluctuation symmetry for a dissipa- 
tion function S has been proposed and studied in various 
contexts. Schematically, that fluctuation relation goes 
somewhat like 

Prob[5 = a] = Prob[S = -a] exper (I) 

where the probability distribution is in the steady state 
and S denotes in general a path-dependent variable 
whose average can often be interpreted as the total 
change (in dimcnsionless units) of entropy over a long 
time interval [0, r]. Depending on the particular real- 
ization, the function S is better called entropy current, 
dissipated work or entropy production. Subject to the 



scale of description S = S v (uj) is a function of the path 
or history oj and of the driving fields ip. We will call S 
the dissipation function. 

Symmetries such as first appeared in 0, Q in the 
context of thermostated and smooth dynamical systems 
where the phase space contraction is identified with en- 
tropy production. A slightly stronger statement is that 
the steady state averages satisfy 

(fe) = {fe- s ) (2) 

where / = / (cj) is a function of the history lu — (u>t)o 
of the system over the time interval [0, r] and 9 is the 
time-reversal, Qlo — (uj T ^ t )o, so that f9(ui) — f(0u). 

While that symmetry (J2J is very generally valid and 
while the symmetry-identity JQ) is non-perturbative (i.e., 
not restricted to close to equilibrium situations), so far 
its further theoretical consequences have only been made 
visible in the derivation of linear response relations such 
as those of Green and Kubo, see e.g. 0, 0|. The bad 
news is that while relations such as JIJ or are very 
general, their consequences for expansions around equi- 
librium are necessarily limited to first order, the linear 
regime, as we show next. 

Expansion of the fluctuation symmetry for the 
dissipation function. As a central assumption, here as 
in irreversible thermodynamics Q we take the dissipation 
function S as given by 

S = S v (uj) = ^2 Va Ja(u) — If ■ J (3) 

a 

where the J Q 's are the currents describing the displace- 
ment of a quantity of type a (matter, charge, energy,...); 
the (p a 's indicate the various accompanying (driving) 
fields or affinities (ip — is set to correspond to equi- 
librium). This assumption implies that we ignore tem- 
poral boundary terms of the form U (u> T ) — U (u>o), as they 
usually play no further role in the fluctuation theory for 
t t +oo, see however 0. We write (-} v for the steady 
state expectation and (-)o = (-)cq stands for the equi- 
librium expectation. No systematic expansion around 
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equilibrium is at all possible from the relations QJ or 
(J2J. That can be seen as follows. 

If /+ is time-symmetric, f + 9 = f + , then (J2J gives to 
linear order in tp, 

(f + ) v = (f + e~ s ) v => (f + ) v = (/+>„ - (/+ tp ■ J)o (4) 

and indeed, the product J/ + is antisymmetric under 
time-reversal, hence vanishes under equilibrium expec- 
tation: results in the identity = to first or- 
der in tp and J2J gives no first order term around equi- 
librium for time-symmetric observables. If / _ is anti- 
symmetric, f~9 = —f~, then the second order term of 
(f~)tp requires information about the first order terms 
(f~J a )tp but that information is lacking because of the 
arguments above applied to the time-symmetric observ- 
able /+ = f~J a - 

Lagrangian set-up. From the above remarks, it is 
clear what is missing. We need information about the 
generation of time-symmetric observables. That is made 
most visible in a Lagrangian set-up. It amounts to the 
realization of the steady state space-time distribution as 
a Gibbs distribution. That can be done under a wide 
variety of contexts, see 0, 0, 0> 13 E5 > but most 
easily and explicitly for stochastic dynamics. 

Quite generally we can construct the probability distri- 
bution P v on space-time configurations to for a nonequi- 
librium steady state parameterized by some vector tp — 
(tp a ), m terms of the corresponding equilibrium distribu- 
tion Pq: 

P v {u)=exp[-C v (u)]I\ ) {u) (5) 

denoting the Lagrangian action by C v to avoid confu- 
sion with the dissipation function S. One cannot go too 
lightly over the choice of the equilibrium distribution Pq 
but we momentarily choose to ignore these issues. 

Relation J5J teaches that the time-antisymmetric part 
is given by the dissipation function: 

S v (w) = L v {Quj) - L v {uj) (6) 

Indeed, with is equivalent with 

P v (u) = e s *Wp v (0u) (7) 

where we used that in equilibrium, Pq9 = Pq; that re- 
lation implies Q). In other words, the source of time- 
reversal breaking is (essentially) given by the dissipation 
function. That has been confirmed by many physically 
motivated examples and models and was also derived in 
much greater generality 0, Q] . It was also argued there 
how JJJ can be interpreted as a unification of existing 
fluctuation theorems, ranging from the Jarzynski equal- 
ity to the Gallavotti-Cohen theorem for chaotic dy- 
namical systems |2|. 

To go beyond linear order in the nonequilibrium 
response-functions, it appears that we must obtain ex- 
tra information about the time-symmetric part in C^. 



It is important to understand here that there is no rea- 
son to think that the nonequilibrium driving would not 
generate an extra term in C v which is symmetric under 
time-reversal. Adding a nonequilibrium driving changes 
the time-symmetric part in the space-time action, which 
is the reason why the response functions cannot be gen- 
erated by the expression for the dissipation function (or 
entropy production) alone. 

Time-symmetric part. For characterizing nonequi- 
librium one has to realize that time-symmetry breaking is 
not spontaneous and is itself very much linked with, if not 
caused by, breaking of spatial or even internal symme- 
tries. A simple case is steady heat conduction made pos- 
sible by the spatial arrangement of different heat baths 
in contact with the system. Therefore, to resolve further 
the nonequilibrium state, one can exploit also directly 
the (anti-)symmetries associated to the driving fields. It 
is straightforward to introduce the notion of field-reversal 
as tp a — > — tp a , which could e.g. mean to exchange the 
two heat reservoirs in a one-dimensional heat conduction 
problem or to reverse the direction of an externally given 
driving field. We utilize field reversal as our basic trans- 
formation on the driving fields and with © 

S- v — —S v , and hence S,p(6u)) = 5_ v ((jj) (8) 

This is different for the time-symmetric part in the La- 
grangian A v = L v + C v 9, for which in general A v ^ 
±A^ V . For Y v = C- v — C v we have the identity 
A- v — Aip = Y v 9 + Y v and we may derive that A v = A- v 
iff Y v = S v iff £- v = C v 9. 

The reason why the time-symmetric term A v is unseen 
in the linear response to currents is because, to linear 
order, field-reversal can be replaced with time-reversal: 
for small ip, 

(J)_ V ^-(J) V = (J9) V 

This is not generally true beyond order tp, which is for ex- 
ample all-important for the construction of so called rec- 
tifiers, see for instance 0|- Hence, while for small driv- 
ing (small tp, close to equilibrium), field-reversal can be 
implemented by time-reversal, the time-symmetric A v is 
expected to be important further away from equilibrium. 
Consider for example the flow of a viscous fluid through a 
tube under influence of a pressure difference tp — Ap. For 
small tp one observes a laminar flow and hence, accord- 
ing to Poiseuille's equation the flow rate is proportional 
to the pressure drop and one cannot distinguish between 
reversal of time and reversal of the field or pressure dif- 
ference. At higher pressure differences however, or above 
a critical velocity, the laminar pattern breaks up and 
the flow becomes turbulent. Then time-reversal of the 
flow pattern definitely differs from the typical pattern ob- 
tained via field-reversal, and hence C v must have a time- 
symmetric part. Ultimately A v ^ A- v (<*=> £- v ^ C v 9) 
is intimately related to nonlinear response. The field- 
antisymmetric part A^ v — A v — Y v 9+Y v rather than the 
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tame-antisymmetric part (the dissipation function S v ) is 
also the pivotal quantity to observe, when one attempts 
higher order expansions. 

New fluctuation symmetry. Remember that Y v — 
£- v — C v . Applying field reversal we have, by construc- 
tion, 



(f)- v = (fe- Y -) v 
Combining it with time-reversal, 



(9) 



(10) 



where we have used that S v = (Y v — Y ip 9)/2. As a con- 
sequence, for R v = (Y v + Y v 9)/2 we have Prob^i?^ = 
r] =Prob_ ¥ ,[i? ¥ , = r]e r . It often happens that the field 
reversal can be implemented as a map on the histories. 
In that case, there is an involution T on paths u> so that 



C- 



C V Y. The transformation T could for example 



simply spatially mirror all the internal degrees of free- 
dom, cf. the examples below. Then, from IjlOjl and simi- 
lar to ©, 



(11) 



As R = R v satisfies ROT = —R, we then have the 
fluctuation symmetry not only for the dissipation 
function S = {Y v — Y v 9)/2 but we now also get it for 
R = (Y v + Y v 6)/2 = {A_ v - A v )/2: 



Prob[i? 



Prob[i? 



-r a 



(12) 



where the probability distribution is in the steady state 
and the fluctuation observable R denotes the variable 
time-symmetric and field-reversal-antisymmetric part in 
the Lagrangian action over a long time-interval [0, r]. 
This fluctuation symmetry is again generally valid, non 
perturbatively and away from equilibrium. This time 
however, it has implications for the response of time- 
symmetric observables. Indeed, from © or l|l(Jf) . when 

/+ = f + o, 



1 

2lp' 
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)tp=Q - 



dip a 



R v )ip=o 



(13) 



As a result, we have an Onsager reciprocity for the ob- 
servables 



d 

dtp c 



R 



(p\ip=0 



(14) 



in the sense that to first order {V a ) v = (V Q )o + 
M a ^(p 7 with symmetric linear response coefficients 
May = M ja ; this can be seen by taking /+ — V 1 in (|13|l . 
The application of i|13|) obviously leads to a higher order 



expansion of the currents J a around equilibrium when 
taking / + = J a J 1 in (JTSJl. 

Examples. We consider here a classical model of heat 
conduction, 0,0]. At each site i — 1, . . . , N there is an 
oscillator characterized by a scalar position qi and mo- 
mentum pi . The dynamics is Hamiltonian except at the 
boundary {1,N} where the interaction with the reser- 
voirs has the form of Langevin forces as expressed by the 
Ito stochastic differential equations 



dq, 
dpi 

dpi 



N 



Pidt, i = 1, 

-T~(q)dt, i = 2 
dqt 



,N-1 



■ — (q)dt-jp i + J 7 fdW i (t), i=l,N 
OQi V Pi 



The Pi , (3n are the inverse temperatures of the heat baths 
coupled to the boundary sites i — 1,N; dWi(t), dWN(t) 
are mutually independent standard white noise. Appro- 
priate growth and locality conditions on the potential U 
(which need not be homogeneous) allow the existence of 
a corresponding smooth Markov diffusion process with a 
unique stationary distribution. When Pi = p^ = P, 
the Gibbs measure ~ exp[— (3H],H — ^2p 2 /2 + U is 
time-reversible for the process with kinematical time- 
reversal 7r given by the involution 7r/(g,p) = f(q,—p). 
Let uj — ((q(t) , p(t)) , t £ [0,t]) denote the evolution of the 
system in the period [0, r]. The natural definition of time- 
reversal 9 is thus (9uj) t — 7r(cj r _ t ) = (q(r — t), — p(r — £)). 
For determining the action in (0 we take for Pq the equi- 
librium process where the inverse temperatures of both 
heat reservoirs are equal to (3 = (Pi + Pm)/2. There is 
one driving field (p = @± — Pn- To compute the anti- 
symmetric part under time-reversal, see ©, in 16] stan- 
dard stochastic calculus yields S v = (Pi — Pn)J where 
J = J(ui) is the (path-dependent) heat dissipated in the 
reservoir at site i = 1 in time t. The steady state average 
of S v is the heat dissipation and is easily seen to be non- 
negative by integrating (JJJ) over all io and by using the 
Jensen inequality, see further in 0, 0, 0] . But we can 
also compute the symmetric part under time-reversal; it 
equals 







p\{t)dt + 



PI 



P 2 



P 



p 2 N {t)dt 



which involves the time-integral of the (path-dependent) 
kinetic temperatures at the ends of the chain. Field re- 
versal corresponds to exchanging Pi with Pn and hence 

R={Pn- Pi) [ I p\{t)dt - f p 2 N (t)dt] 



satisfies the fluctuation symmetry (|12|) . Furthermore, for 
(THI) . we find 

V= [ p 2 N (t)dt- f p\{t)dt 
Ja Ja 
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from which we get that the left /right difference of time- 
integrated kinetic temperatures satisfies a Green-Kubo 
type relation JT3J, but now for the time-symmetric ob- 
servable V = J (p 2 N — Pi)dt, 

A second class of examples is given by the more general 
Langevin-type equation 

dx(t) = -F(x(t)) dt + V2 dW[t) 

The drift is given by the force —F, that, in nonequilib- 
rium situations, cannot be derived from a potential. To 
parameterize the driving, we write F = Vf7 — <pG where 
U is a potential function and G is an external force that 
possibly still depends on the position x(t) of the parti- 
cle. For driving field <p = the process is reversible in 
the stationary distribution ~ exp [—£/]. The path-space 
measure is formally given, [19j, by 

1 f T 

P ip (uj)=exp-- / dt[(x(t)+F(x(t))) - VF(x(t))] 
4 Jo 

The antisymmetric part under time-reversal in the action 
logP^ indeed gives the dissipated power in terms of a 
stochastic Stratonovich integral of the "force" G(x(t)) 
times the "velocity" x(t): S = ip J G{x(t)) x(t) dt. On 
the other hand, the part in C v that is symmetric under 
time-reversal and is antisymmetric under field-reversal is 
given by 

R = ip [ WU(x(t))G(x(t))dt+^ [ VG(x(t))dt (15) 
Jo 2 Jo 

It satisfies the symmetry 1)12(1. When there is a constant 
external field G, the second term vanishes and the first 
term picks up the total force exerted on the particle on its 
trajectory x(t), t S [0, r] . The same can be repeated when 
x(t) = p(r,t) is a field on some bounded domain, r € V, 
that is subject to dissipation and forcing at its 'left' 
and 'right' boundary. The expression of R is likewise of 
the form f gv drdtU'(p(r, t)) — J g v drdtU'(p(r, t)) where 
the spatial integration is over the respective boundaries 
dgV, d r V of V; the integrand U'(p(r,t)) could for exam- 
ple indicate the local chemical potential as defined in the 
system. In that sense the external force couples with the 
local density of the particles. In the case of the previ- 
ous example (heat conduction) one has the external force 
coupling with the local kinetic energy. 

Conclusions. For an irreversible thermodynamics as 
a linear response theory for time-antisymmetric observ- 
ablcs like fluxes or currents, mostly only the entropy pro- 
duction and its fluctuations are needed. Beyond linear 
order, noncquilibrium physics must come to terms also 
with the fluctuations of time-symmetric observables. The 
symmetry that we know is present in the fluctuations of 



the dissipation function is also found to be valid for a 
time-symmetric observable (R). That enables to con- 
tinue the perturbation expansion around equilibrium be- 
yond linear order in the driving fields. As a consequence 
a Green-Kubo type relation and corresponding Onsager 
reciprocity are obtained for time-symmetric functions 
(V). While some examples were shown to give specific 
physical meanings to our R and V, still, noncquilibrium 
effects beyond linear order reside, at least for theory, in a 
vast terra incognita and it will be mainly through exper- 
imental work that the usefulness of our new and general 
relations will be checked. 
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